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基于 微 卫 星 标 记 的 中 国 东 北 地 区 灰 飞 虱 
遗传 变异 及 种 群 遗 传 结构 分 析 
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摘要 :【 目 的 ] 本 研究 旨 在 明确 我 国 东北 地 区 灰 飞 她 Laodelphax striatellus 种 群 遗 传 变异 及 种 群 遗 传 
结构 , 阅 明 种 群 间 遗传 分 化 与 基因 流 。【 方 法 ] 利 用 9 对 微 卫 星 引物 对 采 自 我 国 东北 地 区 15 个 地 
理 种 群 的 375 头 灰 飞 虱 样品 进行 测序 与 分 析 ; 利 用 GeneAlex6.51, GENEPOP4.0.9 和 STRUCTURE 
2.3.4 等 软件 分 析 友 飞 虱 地 理 种 群 间 的 遗传 多 样 性 遗传 分 化 .基因 流 及 种 群 遗传 结构 。[【 结果 】 在 
所 分 析 的 375 头 灰 飞 必 个体 中 ,各 位 点 平均 有 效 等 位 基因 数 Na =6.898; 总 体 上 , 灰 飞 必 不 同 地 理 
种 群 遗 传 多 样 性 较 高 (平均 观测 杂 合 度 Ho =0.548 ;平均 期 望 杂 合 度 He =0.582) ,各 种 群 间 基因 流 
较 低 (Nm =0.660)。UPGMA 聚 类 树 、PCoA 及 STRCTURE 分 析 结 果 表 明 , 东 北 地 区 灰 飞 恒 种 群 分 
为 两 分 支 :吉林 (于 ) 和 沈阳 (SY2012,SY2013 和 SY2014 ) 种 群 聚 为 一 支 ; 其 余 种 群 聚 为 一 支 。 
AMOVA 分 析 结 果 表 明 , 灰 飞 虱 遗传 变异 主要 来 自 种 群 内 (87% ) ,种 群 间 变异 水 平 较 低 (13% ) 。 
【结论 ] 中 国 东北 地 区 灰 飞 虱 遗传 多 样 性 较 高 ,不 同 地 理 种 群 存在 一 定 程度 的 遗传 分 化 , 且 基 因 交 
流 较 低 ,存在 一 定 的 种 群 遗 传 结 构 。 
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中 图 分 类 号 : Q968.1 ”文献 标识 码 : A ”文章 编号 : 0454-6296(2020)01-0073-12 

Genetic variation and population genetic structure of the small brown 
planthopper, Laodelphax striatellus ( Hemiptera: Delphacidae ), in 


Northeast China based on microsatellite markers 

JIANG Shan, WANG Xing-Ya, WANG Xiao-Qi™ (College of Plant Protection, Shenyang Agricultural 
University, Shenyang 110161, China) 

Abstract: [Aim)] The purpose of this study is to clarify the population genetic variation, population 
genetic structure, population genetic differentiation and gene flow of the small brown planthopper, 
Laodelphax striatellus, in Northeast China. [Methods) Nine pairs of microsatellite primers were used to 
sequence and analyze 375 samples of L. striatellus collected from 15 geographical populations in Northeast 
China, and the genetic diversity, genetic differentiation, gene flow and population genetic structure of 
these geographical populations were analyzed using GeneAlex 6.51, GENEPOP 4.0.9 and STRUCTURE 
2.3.4 based on the nine microsatellite (SSR) molecular markers. 【 Results 】 The average number of 
effective alleles per locus was 6. 898 among the 375 individuals of L. striatellus assayed. In general, the 
genetic diversity of different geographical populations of L. striatellus was high ( mean observed 
heterozygosity Ho = 0. 548; mean expected heterozygosity He = 0. 582), and the gene flow between 
different populations was low (Nm = 0. 660). The results of UPGMA clustering tree, PCoA and 
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STRCTURE analysis showed that L. striatellus populations in Northeast China were divided into two 
groups: populations from Jilin (JL) and Shenyang (SY2012, SY2013 and SY2014) were clustered into 
one group, and the remaining populations into another group. The AMOVA analysis showed that the 


genetic variation of L. striatellus was mainly from within the population (87% ), and the inter-population 


variation level was low (13% ). [Conclusion)] The genetic diversity of L. striatellus in Northeast China is 


relatively high, different geographical populations have a certain degree of genetic differentiation, the 


gene exchange is low and a certain population genetic structure exists. 


Key words: Laodelphax striatellus; microsatellite; genetic diversity; genetic structure; gene flow 





灰 飞 恒 Laodelphax striatellus 隶属 于 半 怒 日 
( Hemiptera) 飞 乱 科 (Delphacidae ) ,是 我 要 的 农 
业 害 忠 之 一 ,可 人 危害 水 称 、 麦 类 、 玉 米 等 多 种 不 本 科 
植物 , 常 可 造成 严重 的 经 济 损失 ( 顾 伯 良 等 , 2005 ) 。 
该 种 害虫 广泛 分 布 于 东亚 、 东 南亚 、 欧 洲 和 北美 等 地 
区 , 尤 以 东亚 温带 地 区 发 生 危 害 最 为 严重 ,在 我 国 水 
稳 主 产 区 普遍 分 布 (刘向东 等 ,2006)。20 世纪 60 
年 代 , 该 种 害虫 在 江苏 大 面积 暴发 ,而 后 80 年 代 暴 
发 趋势 蔓延 到 山东 。 近 年 来 , 随 着 免 耕 、 稳 套 麦 及 麦 
套 稻 等 轻型 栽培 技术 以 及 感 病 高 产品 种 的 推广 应 用 
( 印 德 良 等 , 2005 ) ,同时 受 冬季 平均 气温 的 升 高 以 
及 生态 环境 的 恶化 等 诸多 因素 的 影响 ( 吴 雪 分 等 ， 
2005 ) , 灰 飞 翅 的 发 生 量 逐 年 增加 。 特 别 是 近年 来 ， 
由 于 长 期 大 量 使 用 化 学 杀 虫 剂 ,使 得 灰 飞 乱 对 一 些 
药剂 产生 了 严重 抗 药性 ,给 该 种 害虫 的 防治 带 来 更 
大 困难 ( 孙 志 广 等 , 2018 ) 。 

微 卫 星 标记 (microsatellite marker) ,又 被 称 为 简 
单 序列 重复 (simple sequence repeats，SSR ) ,是 种 群 
遗传 学 和 分 子 生态 学 等 研究 领域 的 重要 研究 手段 ， 
具有 高 多 态 性 .自然 选择 压力 小 、 共 显 性 遗传 .实验 
重复 性 高 .对 DNA 模板 要 求 低 以 及 突变 速率 高 等 优 
点 而 广泛 应 用 于 遗传 图 谱 构 建 、 克 隆 、 农 作物 育种 、 
分 子 标记 筛选 .分子 谱系 地 理 研 究 等 领域 ( Jarne and 
Lagoda, 1996; Tang et al., 2015; Duan et al., 
2017)。 此 外 ,研究 物种 种 群 遗传 结构 是 害虫 综合 
治理 的 重要 内 容 , 有 助 于 探讨 其 种 群 间 遗 传 进化 关 
系 ,揭示 其 起 源 、 扩 散 路 线 及 趋势 ,进而 为 明确 该 种 
害虫 的 发 生动 态 与 成 灾 适 应 机 制 提供 理论 基础 。 

近年 来 , 随 着 生物 技术 的 进步 , 微 卫 星 的 开发 和 
分 型 技术 越 来 越 普遍 。 目 前 ,我 国 对 于 灰 飞 恒 的 种 
群 遗传 多 样 性 种群 扩 散 及 历史 动态 等 方面 研究 已 
有 部 分 报道 。Hoshizaki (1997 ) 利用 同 工 酶 技术 对 
日 本 和 我 国 台湾 地 区 不 同 地 理 种 群 的 灰 飞 乱 进 行 遗 
传 多 态 性 分 析 ,发 现 不 同 地 区 种 群 灰 飞 乱 存在 地 理 
种 群 差 异 ; 万 由 衷 等 (2001 ) 对 不 同 地 理 种 群 的 灰 飞 































































































异 进 行 RAPD 分 析 ,将 其 分 为 两 个 复 , 禾 1 包括 北 
京 、 上 海 、 福 建 海南 .辽宁 、 四 川 及 云南 种 群 ; 自 2 包 
括 宁 夏 和 日 本 种 群 。 并 且 不 同 地 区 灰 飞 乱 种 群 的 遗 
传 距离 与 其 地 理 距 离 不 呈 相 关 性 。 对 云南 省 和 江苏 
省 中 7 个 灰 飞 乱 的 地 理 种 群 研究 表明 , 灰 飞 乱 rDNA 
的 基因 间隔 区 (intergenic spacer，IGS ) 在 种 内 存在 
多 态 性 ,7 个 灰 飞 各 种 群 未 形成 明显 的 地 理 种 群 分 
化 ( 座 富 荣 , 2004) 。 利 用 642 bp 的 线粒体 基因 组 
片段 对 15 个 地 理 种 群 309 头 灰 飞 乱 的 测序 结果 进 
行 分 析 , 共 确定 16 种 单 倍 型 ,其 中 , 单 倍 型 1 和 单 倍 
型 2 是 最 主要 的 单 倍 型 , 占 整体 的 87.7% ,不 同 地 
理 种 群 灰 飞 避 遗传 分 化 程度 较 小 (Zhang et al.， 
2013 ) 。 利 用 线粒体 基因 (CO7 和 CON) 和 13 个 SSR 
位 点 对 中 国 22 个 灰 飞 乱 地 理 种 群 进行 种 群 遗 传 结 
构 及 种 群 遗传 多 样 性 的 研究 发 现 ,整个 气候 区 的 微 
卫星 多 样 性 水 平 大 致 相似 (Sun et al., 2015 ) 。 

本 研究 利用 微 卫 星 标记 对 中 国 东 北 地 区 灰 飞 恒 
不 同 地 理 种 群 的 遗传 多 样 性 和 遗传 结构 进行 研究 ， 
以 期 分 析 中 国 东 北 地 区 灰 飞 虱 种 群 间 的 遗传 分 化 和 
基因 流 , 为 制定 灰 飞 恒 综合 防治 措施 在 分 子 生物 学 
方面 提供 借鉴 。 










































































1 材料 与 方法 


1.1 供 试 虫 源 

2012 - 2016 年 ,实地 采集 我 国 东 北 13 个 市 县 
的 375 涉 灰 飞 恒 成 虫 样品 ( 表 1) ,其 中 包括 分 别 在 
2012 -2014 年 采集 沈阳 3 个 种 群 。 采 集 时 每 头 灰 
飞 乔 成虫 保持 3 m 以 上 距离 ,避免 采集 到 来 自 同一 
个 父母 本 后 代 个 体 。 将 收集 到 样品 浸泡 在 -20%C， 
95% 乙醇 中 ,保存 于 沈阳 农业 大 学 植物 保护 学 院 。 
1.2 灰 飞 虱 基因 组 DNA 提取 

将 单 头 试 虫 置 于 1.5 mL 离心 管 中 , 加 入 蒸馏 水 
多 次 洗涤 后 晾 干 ,使 用 天 根 生 化 (北京 ) 有 限 公 司 的 
血液 /细胞 /组 织 基 因 组 DNA 提 取 试 剂 盒 提取 灰 飞 
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表 1 中 国 东北 地 区 灰 飞 虱 采集 信息 及 供 试 个 体 数量 
Table 1 Collecting data of Laodelphax striatellus in Northeast China 
种 群 代码 ; 样本 数 虫 态 采集 时 间 ( 年 -月 -日 ) 
采集 地 点 经 纬度 
Population Sample Developmental Collecting date 
Collecting locality Latitude and longitude 
code size stage (year-month-day ) 
辽宁 凤 城 40°23'48.08”N, 
FC 30 成 虫 Adult 2015-05-13 
Fengcheng，Liaoning 124°04'18.05"E 
辽宁 盘锦 41°13'54.61”N, 
PJ 19 成 虫 Adult 2015-04-22 
Panjin, Liaoning 121°59'49.25"E 
辽宁 大 石 桥 40°38'54.41”N, 
DSQ 31 成 虫 Adult 2015-05-13 
Dashiqiao, Liaoning 122°31'13.04"E 
辽宁 沈阳 沈河 41°49'26. 44"N， 
SY2012 22 成 虫 Adult 2012-04-27 
Shenhe, Shenyang, Liaoning 123°32'44.11"E 
辽宁 沈阳 沈河 41°49'26.44"N， , 
SY2013 10 成 虫 Adult 2013-04-27 
Shenhe, Shenyang, Liaoning 123°32'44.12"E 
辽宁 沈阳 沈河 41°49'26.44"N， > 
SY2014 13 成 虫 Adult 2014-04-27 
Shenhe, Shenyang, Liaoning 123°32'44.13"E 
了 宁 台 安 41°23'14.44”N， 
AS 本 全 全 27 成 虫 Adult 2015-05-13 
Tai'an，Liaoning 122°25'50. 65"E 
辽宁 普兰 店 39°24'3.57"N, 
PLD 站 16 成 虫 Adult 2015-05-13 
Pulandian, Liaoning 121°57'39.92"E 
辽宁 开 原 42°32'46.71”N，, 
KY 30 成 虫 Adult 2014-09-03 
Kaiyuan, Liaoning 124°02'17.76"E 
辽宁 新 民 41°59'6.69"N,， 
HT 29 成 虫 Adult 2015-05-13 
Xinmin, Liaoning 122°50'12.21"E 
黑龙 江 牲 丹 Y 44°35'37. 86"N， 
MDJ 30 成 虫 Adult 2016-08-21 
Mudanjiang，Heilongjiang 129°37'21.60"E 
黑龙 江 安 达 46"25'27.05"N， 
AD 43 成 虫 Adult 2016-08-22 
Anda，Heilongjiang 125°20'15. 39"E 
吉林 白城 45°37'8.95"N, 
BC 15 成 虫 Adult 2016-09-05 
Baicheng, Jilin 122°51'45.60"E 
吉林 吉林 43°50'44.09"N， 
JL 30 成 虫 Adult 2016-08-20 
Jilin，Jilin 126°33'17.94"E 
吉林 松原 45°08'27.23"N, 
SYS 30 成 虫 Adult 2016-08-22 
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恒 总 DNA。 用 NanoDrop ND-2000 超 微量 核酸 和 蛋白 
测定 仪 检测 提取 的 灰 飞 乱 基因 组 DNA 浓度 以 及 吸 
光度 值 ,将 检测 合格 的 DNA 样品 置 于 -20%C 冰箱 内 
保存 。 
1.3 微 卫 星 位 点 的 筛选 

根据 孙 研 涛 (2012) 和 Sun 等 (2015 ) 的 16 个 灰 
《和 翅 微 卫星 标记 进行 筛选 , 微 卫 星 引 物 由 上 海 生 工 
生物 工程 技术 服务 有 限 公 司 合成 。 经 过 筛选 ,从 中 
选取 具有 多 态 性 高 ,易于 扩 增 的 9 个 微 卫 星 位 点 , 引 
物 信息 列 于 表 2。 

PCR 反应 体系 (20 pL): ddH,0 16.3 pL, 2.5 








mmol/L dNTPs 0.2 pL, 10 x EasyTaq Buffer 2 AL， 
上 下 游 引物 (10 pmol/L) 各 0.4 nL, EasyTaq DNA 
Plymerase 0.2 pL, DNA 模板 0.5 LLo 反应 条 件 : 
94°C 4 min; 94°C 30 s, 53 ~ 58°C 30 s, 72°C 30 s， 
共 35 个 循环 ; 72% 延伸 5 min, 4% 保存 。PCR 扩 增 
产物 通过 1.0% 琼脂 糖 凝 胶 电 泳 检测 扩 增 情况 后 ， 
将 目标 产物 送 往 上 海 生 工 生物 技术 有 限 公 司 进行 毛 
i 最 后 用 GeneMarker 软件 
行 等 位 基因 片段 大 小 的 读 取 ,并 人 工 进 行 位 点 数 
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表 2， 微 卫星 引物 信息 


Table 2 Information of microsatellite primers 








Annealing temperature 


退火 温度 (%C ) GenBank 登录 号 扩 增 长 度 (bp) 


Amplified length 


GenBank accession number 





位 点 引物 序列 (5' 31) 重复 单元 

Locus Primer sequence Repeat motif 
AGAGAGAGAGAGAGACACAC 

LS1 (AC)s, (AG)se 
GAAAAAGCACTTGCCACATT 
TCTCTCTCTCTCTCACACAC 

LS3 (AC)s 
GCGGTCGCTAATACACTCC 
AGAGAGAGAGAGAGACACAC 

LS6 (AC)7 
TAATACAGGGTGCGTCGTTAT 
AGAGAGAGAGAGAGACACAC 

LS7 (AC)s 
CTACCATCCATCGGAATGG 
TCTCTCTCTCTCTCACACAC 

LS8 《AC) mn 


AACTCATTTCATAGCCCCAAC 


GCACAGTGAAATTGCGTCAT 
LS13 (GT), 
TGTGAGTGAATTCGGGTGAA 


TGCACTGACACCACAGTCAA 
LS12 (GGC)s 
AAGGTGTCGTCGTACCTGCT 


TCCAGCTACGGACAACAACA 
LS15 (ACA)6 
AAGGTGTCGTCACCCATCAT 


CCACTGGGATAAATCGATGG 
LS16 (GCA); 
AGGATGAGCGGACAAAGAGA 


1.4 数据 分 析 

本 研究 利用 GENEPOP 4. 0. 9 (Raymond and 
Rousset，1995 ) 进行 灰 飞 乱 各 个 位 点 和 种 群 哈 迪 - 温 
伯 格 平衡 检测 及 连锁 不 平衡 检测 。 利 用 GenAlex 
6. 51( Peakall and Smouse,，2006 ) 软件 计算 种 群 在 各 
位 点 上 的 等 位 基因 数 ( Ne) 、 观 测 杂 合 度 (Ho) 、 期 望 
杂 合 度 (He) 、 非 偏差 期 户 杂 合 度 (wuHe)、Shannon 信 
息 指 数 (1) 、F-statistics 统计 和 基因 流 (Nm) ,以 及 位 
点 在 各 种 群 中 的 等 位 基因 数 、 观 测 杂 合 度 ,期望 杂 合 
度 以 及 Nei 氏 期 望 杂 合 度 。 并 根据 遗传 距离 矩阵 进 
行 PCoA (principal coordinates analysis ) 分 析 。 使 用 
MEGA7( Kumar et al., 2016) 基 于 Nei 氏 遗 传 距离 运 
用 UPGMA 聚 类 法 构建 系统 发 生 树 。 使 用 
Arlequin3.5.2.2( Excoffier et ol.，2005 ) 软件 分 析 种 
群 间 遗 传 分 化 和 基因 流水 平 . 用 STRUCTURE 2.3.4 
(Pritchard et al., 2000) 根 据 贝 叶 斯 聚 类 的 方法 分 析 
种 群 遗 传 结构 ,重复 10 000 次 。 使 用 Evanno 的 AK 
方法 (Evanno et al.， 2005 ) 确定 最 佳 的 K 值 后 使 用 
CLUMPP1. 1. 2 软件 (Jakobsson and Rosenberg， 
2007 ) 对 结果 进 比 对 和 整合 并 使 用 Excel 2016 作 图 。 
使 用 CeneAlex 6. 51 软件 进行 分 子 方差 分 析 
(AMOVA) ,将 15 个 灰 飞 乱 地理 种 群 作为 整体 ,按照 
贝 叶 斯 聚 类 结果 分 别 进行 分 组 ,然后 分 别 在 这 两 种 
























































53 JN835260 97 -177 
53 JN835262 201 -259 
53 JN835265 126 -147 
53 JN835266 91 - 123 
55 JN835267 84 -142 
57 SRR035224. 4952 134 -166 
58 SRR035224. 123990 181 -199 
57 SRR035224. 153340 92 -141 
56 SRR035224. 161832 96 -141 


分 组 方式 下 进行 AMOVA 分 析 。 
2 结果 


2.1 灰 飞 恤 种 群 微 卫星 位 点 的 检测 
通过 荧光 引物 优化 和 筛选 ,确定 了 9 对 扩 增 效 
果 较 好 且 多 态 性 高 的 微 卫 星 位 点 ( 表 2) 。 并 按理 论 
片段 大 小 将 9 对 荧光 标记 分 为 3 组 。 使 用 
GENEPOP 4.0.9 对 各 位 点 进行 哈 迪 - 温 伯 格 平衡 检 
验 , 其 结果 经 过 Bonferroni 校正 ,大 部 分 位 点 未 偏离 
哈 迪 - 温 伯 格 平衡 (P < 0.00024) , 且 偏 离 哈 迪 - 温 伯 
格 平衡 的 种 群 基本 都 表现 出 杂 合 子 缺 失 情 况 ( 表 
3) 。 对 微 卫星 之 间 的 连锁 不 平衡 性 测验 发 现 , 各 位 
点 间 不 存在 显著 连锁 现象 (P <0.05)。 

9 个 位 点 无 效 等 位 基因 频率 为 0.021 ~0. 127， 
其 中 位 点 LS1 的 无 效 等 位 基因 频率 最 高 为 0. 127， 
所 有 位 点 的 无 效 等 位 基因 频率 均 小 于 0.200。 由 中 
国 东 北 地 区 15 个 灰 飞 乱 地 理 种 群 在 9 个 微 卫 星 位 
点 上 的 遗传 多 样 性 参数 值 ( 表 4) 可知 ,每 个 位 点 的 
平均 等 位 基因 数 为 20 ,平均 观测 杂 合 度 为 0. 566 , 平 
均 期 望 杂 合 度 为 0.814 ,平均 Nei 氏 期 望 杂 合 度 为 
0.813。 总 体 上 , 灰 飞 乔 各 位 点 等 位 基因 数 较 高 , 杂 
合 度 水 平 较 高 且 近 交 系数 低 ,表明 遗传 多 样 性 较 高 。 
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表 3 9 个 微 卫 星 位 点 在 中 国 东北 地 区 各 灰 飞 虱 种 群 中 的 哈 迪 - 温 伯 格 平衡 的 检测 

Table 3 Hardy-Weinberg equilibrium test in each population of Laodelphax striatellus in Northeast China at nine SSR Locus 
种 群 代码 位 点 Locus 
Population code LS1 LS3 LS6 LS7 LS8 LS13 LS12 LS15 LS16 
FC 0.002 0.001 0.000” 0.000” 0.004 0.855 0.001 0.001 0. 122 
PJ 0.019 0.000” 0.001 0.227 0.017 0.031 0.035 0.000” 0.000” 
DSQ 0.028 0.000” 0.001 0.000” 0.024 0.163 0.991 0.001 0.670 
SY2012 0.242 0.929 0.000” 0.001 0.230 0.000” 0.042 0. 129 0.000” 
SY2013 0.000” 1.000 0.000” 1.000 0.084 0.000” 1.000 0.000”* 0.007 
SY2014 0.004 0.845 0.000” 0.729 0. 004 0.000” 1. 000 0. 122 0.002 
AS 0.006 0.264 0.367 0.499 0.217 0.119 0.100 0.838 0.224 
PLD 0.014 0.249 0.245 0.753 0.166 0.518 0.423 0.293 0.628 
KY 1.000 0.000” 0.025 0.818 0.017 1.000 0.264 0.158 1.000 
HT 0.000” 0. 134 0.027 0.988 0.001 0.612 0.000” 0.855 0.095 
MDJ 0.006 0.000” 0.000” 0.012 0.088 0.444 0.310 0.236 0.000” 
AD 0.001 0.594 0.030 0.005 0.000” 0.348 0.073 0.000” 0.000” 
BC 0.031 1.000 1.000 0.554 0.061 0.093 0.830 0.434 0.570 
]L 0.011 0.851 0.000” 0.018 0.031 0.000” 0.017 0. 198 0.017 
SYS 0.000” 0.032 0.562 0.219 0.000” 0.854 0.022 0.001 0.461 

种 群 代码 见 表 1 ;下 同 。For population code，see Table 1. The same below， 星 号 表示 经 Bonferroni 校正 后 与 哈 迪 温 伯 格 平衡 存在 显著 偏差 (已 < 


0. 00024 ) 。 








Asterisk indicates significant deviation from Hardy-Weinberg equilibrium after Bonferroni correction (P<0.00024). 


表 4 中 国 东北 地 区 15 个 灰 飞 虱 地 理 种 群 中 9 个 微 卫星 位 点 的 遗传 多 样 性 


Table 4 Genetic diversity of nine microsatellite loci in 15 populations of Laodelphax striatellus in Northeast China 


有 效 等 位 








Shannon 














、 ”等 位 基因 数 上 观测 杂 合 度 期望 杂 合 度 Nei 氏 期 望 无 效 等 位 
Rs 样本 量 基因 数 信息 指数 x a 
位 点 Number of l Observed Expected 合 度 基因 频率 
Sample Number of effective & Shannon-Weiner a . 
Locus alleles : heterozygosity heterozygosity Nei’s expected Null allele 
size alleles index . 
Na Ho He heterozygosity frequency 
Ne I 
LS1 654 34 3.919 2.130 0.312 0.746 0.745 0.127 
LS3 646 45 14.303 3.071 0.678 0.932 0.930 0.042 
LS6 548 13 4.799 1.888 0.460 0.793 0.792 0.091 
LS7 706 17 6.932 2.204 0.629 0.857 0.856 0.042 
LS8 690 33 11.563 2:753 0.641 0.915 0.914 0.084 
LS13 596 8 2.399 1.198 0.456 0.584 0.583 0.037 
LS12 714 9 5.250 1.823 0.487 0.811 0.810 0.101 
LS15 728 12 6.459 2.067 0.714 0. 846 0. 845 0.031 
LS16 744 16 6.452 2. 173 0.712 0. 846 0. 845 0.021 
平均 值 Mean 670 20 6.898 2. 145 0.566 0.814 0.813 0.064 
此 外 ( 表 5) ,各 个 灰 飞 乱 种 群 平 均等 位 基因 数 为 下 。9 个 微 卫星 位 点 在 15 个 灰 飞 乱 地 理 种 群 内 近 
5.511 ,平均 有 效 等 位 基因 数 为 3. 253 ,平均 观测 杂交 系数 (mrx) 在 -0.223 ~0.367 之 间 , 群体 总 近 交 
合 度 为 0. 548 ,平均 期 望 杂 合 度 为 0.582。Shannon ”系数 (了) 范围 为 0.098 ~0.618。 种 群 间 分 化 系数 


言 息 指数 在 0. 614 ~ 1. 679 间 。 其 中 , 开 原 (KY ) 种 
群 遗 传 多 样 性 最 低 ,普兰 店 (PLD ) 种 群 遗 传 多样 性 
则 最 高 。 

2.2 灰 飞 虱 种 群 遗传 分 化 与 基因 流 

通过 每 个 位 点 的 固定 指数 已 ,， 巴 和 开 w 检 验 种 
的 遗传 分 化 ,各 个 位 点 的 F-statistics 分 析 结 果 如 









































(Fw) 范 围 为 0. 178 ~ 0.462。 根 据 基因 流 的 大 小 ， 
， LS8 的 基因 交流 程度 相对 较 大 (1 < Nm <4)， 
其 余 位 点 间 基 因 交 流程 度 较 小 (Nm <1) ,总体 上 ， 
各 位 点 间 基 因 流 大 小 为 0.660。 

灰 飞 异 不 同 地 理 种 群 间 遗传 分 化 (F's) 及 基因 
流 (Nm ) 研究 结果 表明 , 各 地 理 种 群 间 的 遗传 分 化 指 
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表 5 中 国 东北 地 区 15 个 灰 飞 虱 地 理 种 群 的 遗传 多 样 性 
TableS Genetic diversity among 15 populations of Laodelphax striatellus in Northeast China 
有 效 等 位 
等 位 基因 数 ”基因 数 Damn。 观测 杂 合 度 。 期望 杂 合 度 一 二 生 。 。 转 定 指数 
种 群 代码 Number of Number of 信息 和 数 Observed Expected . 杂 合 度 Fixation 
Population code alleles effective ss heterozygosity heterozygosity 人 ee index 
人 ae index H, 六 heterozygosity 关 
i 了 uHe 

FC 7.111 +0.841 4.199 +0.614 1.516 +0.185 0.573 +0.082 0.682+0.073 0.693 +0.0734 0.168 +0.096 
PJ 6.778 +0.812 3.570 +0.358 1.464 +0.114 0.550 +0.090 0.688 +0.044 0.714+0.046 0.204+0.116 
DSQ 7.778 +1.382 4.785 +0.988 1.569 +0.211 0.635 +0.085 0.696 +0.075 0.708 +0.076 0.100+0.067 
SY2012 4.667 +1.202 2.580 +0.679 1.025 +0.246 0.582+0.128 0.496+0.102 0.516 +0.108 -0.213 +0.154 
SY2013 3.111 +0.790 1.977 +0.489 0.705 +0.217 0.437 +0.135 0.358 +0.103 0.394 +0.114 -0.211 +0.118 
SY2014 4.222 +1.051 2.815 +0.865 0.972 +0.246 0.469+0.122 0.473+0.104 0.505 +0.113 -0.063 +0.186 
AS 4.556 +0.530 2.956 +0.382 1.160 +0.120 0.584+0.067 0.616 +0.0463 0.627 +0.047 0.067 +0.060 
PLD 7.890 +1.060 4.838 +0.695 1.679 +0.149 0.719+0.066 0.749 +0.043 0.774 +0.044 0.055 +0.057 
KY 2.890 +0.389 1.638 +0.140 0.614 +0.093 0.374 +0.071 0.354 +0.055 0.360 +0.056 -0.045 +0.085 
HT 6.333 +0.553 3.849 +0.471 1.440 +0.130 0.614+0.098 0.692 +0.053 0.705 +0.054 0.134+0.116 
MDJ 5.890 +0.696 3.531 +0.507 1.352 +0.151 0.611 +0.078 0.658 +0.054 0.670 +0.055 0.061 +0.010 
AD 7.444 +0.899 4.313 +0.605 1.551 +0.155 0.546 +0.061 0.708 +0.059 0.718 +0.059 0.218 +0.065 
BC 4.000 +0.408 2.406 +0.224 0.993 +0.091 0.560 +0.038 0.557+0.040 0.577 +0.041 -0.018 +0.051 
J]L 4.000 +1.014 1.823 +0.418 0.665 +0.225 0.333 +0.128 0.319+0.110 0.328 +0.112 0.241 +0.162 
SYS 6.000 +0.553 3.511 +0.442 1.374+0.114 0.630+0.064 0.679 +0.038 0.690 +0.039 0.074 +0.079 
We 5.511 +0.812 3.253 +0.525 1.205 +0.163 0.548 +0.088 0.582 +0.067 0.599 +0.069 0.051 +0. 101 


数 为 0.012 ~ 0.550。 其 中 ,沈阳 (SY2012) 和 沈阳 
(SY2014) 种 群 间 遗传 分 化 最 小 (Fs =0.012) ,吉林 
(JL) 和 开 原 (KY ) 种群 间 遗传 分 化 最 大 (Fs = 
0.550)。 各 地 理 种 群 间 基因 流 Nm 为 0. 204 ~ 
20. 275。 其 中 ,沈阳 (SY2012) 和 沈阳 (SY2014 ) 种 
群 间 基因 流 最 大 (Nm =20. 075 ) ,吉林 ( 呈 ) 和 开 原 


表 6 














CKY ) 种 群 间 基 因 流 最 小 (Nm =0. 204)。 此 外 ,吉林 
(L) 种 群 与 其 他 14 个 地 理 种 群 的 基因 流 (Nm) 均 
小 于 1, 遗传 分 化 系数 (Fs ) 均 大 于 0.25。 因 此 , 吉 
林 种 群 与 其 他 种 群 存在 较 大 的 遗传 分 化 ,其 他 大 部 
分 种 群 间 遗 传 分 化 不 明显 ( 表 6) 。 





























中 国 东北 地 区 15 个 灰 飞 虱 地 理 种 群 间 成 对 固定 系数 Psr( 下 三 角 ) 及 基因 流 Nm( 上 三 角 ) 分 析 


Table 6 Pairwise Fs ( below diagonal) and gene flow Nm ( above diagonal) among 1S populations of 
Laodelphax striatellus in Northeast China 





a oie FC PJ DSQ SY2012 SY2013 SY2014 AS PLD KY HT MDJ AD BC JL SYS 
FC 1.014 5.671 0.793 0.560 0.634 1.098 1.593 1.593 0.979 1.197 1.339 0.740 0.334 1.352 
PJ 0. 198 0.967 0.992 0.696 0.806 0.856 1.757 0.356 0.991 1.268 1.339 0.740 0.332 1.029 
DSQ 0.042 0.205 0.734 0.536 0.611 1.023 1.584 0.377 1.051 1.330 1.236 0.741 0.371 1.287 
SY2012 0.240 0.201 0.254 2.940 20.273 0.668 1.399 0.395 0.976 0.732 1.029 0.862 0.358 0.947 
SY2013 0.309 0.264 0.318 0.078 10.051 0.438 0.784 0.235 0.604 0.504 0.636 0.496 0.252 0.553 
SY2014 0.283 0.237 0.290 0.012 0.024 0.530 1.097 0.294 0.813 0.634 0.757 0.563 0.339 0.660 
AS 0.185 0.226 0.196 0.272 0.363 0.321 2.036 0.425 1.331 1.317 1.494 0.903 0.315 2.041 
PLD 0.136 0.125 0.136 0.152 0.242 0.186 0.109 0.735 3.326 2.480 2.480 1.157 0.535 3.469 
KY 0.402 0.412 0.399 0.387 0.516 0.460 0.370 0.254 0.536 0.504 0.477 0.257 0.204 0.487 
HT 0.203 0.201 0.192 0.204 0.293 0.235 0.158 0.070 0.318 1.347 1.332 1.009 0.537 1.298 
MDJ 0.173 0.165 0.158 0.255 0.331 0.283 0.160 0.092 0.332 0.157 1.723 0.713 0.348 1.638 
AD 0.157 0.141 0.168 0.195 0.282 0.248 0.143 0.092 0.344 0.158 0.127 0.904 0.461 2.352 
BC 0.253 0.255 0.252 0.225 0.335 0.308 0.217 0.178 0.493 0.199 0.260 0.217 0.316 1.041 
JL 0.428 0.430 0.402 0.411 0.498 0.424 0.442 0.319 0.550 0.318 0.418 0.352*0.442”* 0.468 
SYS 0.156 0.195 0.163 0.209 0.311 0.275 0.109 0.067 0.339 0.162 0.132 0.096 0.194 0.348 


* P<0.05( 多 重 检验 Multiple test). 
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区 灰 飞 乔 遗传 变异 及 种 群 遗 传 结构 分 析 





研究 结果 表明 ,不 同 地 理 种 群 灰 飞 乱 的 遗传 距 
离 在 0. 135 ~ 2. 198 之 间 , 凤 城 (FC) 和 大 石 桥 
(DSQ) 种 群 间 的 遗传 距离 最 小 (0. 135) ,吉林 (JL) 
种 群 和 凤 城 (FC) 种 群 间 遗 传 距离 最 大 (2. 198)( 表 


表 7 














7)。 此 外 ,不 同 地 理 种 群 灰 飞 乱 的 遗传 相似 度 为 
0. 111 ~0. 874, 风 城 (FC) 和 吉林 (于 ) 种 群 间 遗 传 相 
似 度 最 小 (0.111) , 凤 城 (FC) 和 大 石 桥 (DSQ ) 种 群 
遗传 相似 度 最 大 (0.874) 。 











中 国 东北 地 区 15 个 灰 飞 虱 地 理 种 群 间 的 遗传 距离 (下 三 角 ) 和 遗传 相似 度 ( 上 三 角 ) 


Table 7 Genetic distance ( below diagonal) and genetic similarity ( above diagonal) 


among 15 populations of Laodelphax striatellus in Northeast China 




















种 群 代码 FC PJ DSQ SY2012 SY2013 SY2014 AS PLD KY HT MDJ AD BC JL SYS 
Population code 
FC 0.398 0.874 0.252 0.171 0.176 0.544 0.553 0.264 0.392 0.538 0.527 0.379 0.111 0.567 
PJ 0.920 0.352 0.303 0.231 0.283 0.402 0.551 0.305 0.366 0.541 0.547 0.360 0.123 0.406 
DSQ 0.135 1.044 0.211 0.161 0.158 0.499 0.534 0.253 0.414 0.570 0.474 0.361 0.140 0.533 
SY2012 1.379 1.193 1.554 0.715 0.785 0.221 0.475 0.399 0.357 0.224 0.240 0.399 0.201 0.308 
SY2013 1.769 1.464 1.827 0.336 0.833 0.129 0.345 0.292 0.244 0.135 0.137 0.343 0.129 0.190 
SY2014 1.737 1.263 1.842 0.242 0.183 0.198 0.463 0.340 0.353 0.220 0.213 0.358 0.190 0.264 
AS 0.609 0.910 0.696 1.510 2.048 1.618 0.717 0.446 0.613 0.637 0.640 0.561 0.127 0.750 
PLD 0.593 0.596 0.627 0.745 1.064 0.770 0.332 0.678 0.761 0.728 0.725 0.547 0.342 0.784 
KY 1.334 1.187 1.374 0.919 1.233 1.079 0.807 0.388 0.501 0.494 0.378 0.221 0.283 0.455 
HT 0.937 1.004 0.883 1.030 1.412 1.042 0.490 0.273 0.690 0.577 0.515 0.531 0.384 0.539 
MDJ 0.621 0.615 0.562 1.498 2.006 1.512 0.452 0.318 0.705 0.550 0.649 0.393 0.166 0.660 
AD 0.641 0.603 0.746 1.427 1.986 1.545 0.447 0.322 0.973 0.663 0.432 0.456 0.243 0.728 
BC 0.970 1.023 1.019 0.920 1.070 1.026 0.577 0.603 1.511 0.634 0.935 0.785 0.185 0.559 
JL 2.198 2.095 1.967 1.607 2.045 1.661 2.066 1.072 1.262 0.958 1.794 1.415 1.686 0.266 
SYS 0.567 0.900 0.629 1.177 1.660 1.333 0.288 0.244 0.788 0.618 0.416 0.318 0.582 1.324 

2.3 灰 飞 虱 地 理 种 群 遗传 结构 基于 遗传 距离 矩阵 构建 的 PCoA 结果 显示 ,第 1 

利用 UPGMA 法 对 15 个 灰 飞 乱 地 理 种 群 进行 和 第 2 坐标 轴 上 SSR 变异 贡献 率 分 别 为 43% 和 





况 站 


3 UN 一 


Nei 氏 聚 类 分 析 , 在 相似 系数 为 0. 80 处 15 个 
灰 飞 乱 种 群 被 两 支 ,第 一 支 为 吉林 市 ( 开 ) 种群 ,第 
二 支 为 其 他 种 群 。 而 在 相似 系数 为 0.70 处 ,第 二 支 
种 群 又 被 分 为 沈阳 种 群 (SY2012，SY2013 和 
SY2014) 和 其 他 种 群 (图 1)。 
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中 国 东 北 地 区 15 个 灰 飞 二 地 理 种 群 间 基于 
Nei 氏 遗 传 距离 的 UPGMA 聚 类 图 
UPGMA dendrogram for 15 geographical populations 


0.80 
图 1 














Fig. 1 
of Laodelphax striatellus in Northeast China based on 


Nei’ s genetic distance 














14% ,基于 贝 叶 斯 聚 类 的 中 国 东北 地 区 灰 飞 乱 15 个 
地 理 种 群 也 产生 了 类 似 的 遗传 结构 (图 2) ,其 中 沈 
阳 种 群 (SY2012，SY2013 和 SY2014 ) 和 吉林 种 群 
( 开 ) 倾 向 于 聚 为 一 组 ,而 其 余 种 群 种 群 聚 为 一 组 。 

应 用 STRUCTURE 2.3.4 软件 对 每 个 值 均 运 
算 10 次 , 当 K=2 时 (图 3) ,为 最 佳 分 组 方式 。 聚 类 
结果 显示 , 当 K =2 时 ,吉林 种 群 ( 开 ) 沈阳 种 群 
(SY2012, SY2013 和 SY2014) 和 开 原 种 群 (KY) 聚 
为 一 个 类 群 ,其 余 种 群 聚 为 一 个 类 群 (图 4) 。 
2.4 AMOVA 分 子 方差 分 析 

本 研究 将 中 国 东北 地 区 灰 飞 乱 种 群 作为 整体 进 
行 分 子 方差 分 析 ( 表 8) ,根据 AK 值 将 15 个 地 理 种 
群 分 为 2 个 组 ,综合 PCoA 结果 ,第 1 组 为 沈阳 种 群 
(SY2012, SY2013 和 SY2014) 和 吉林 种 群 (JL) 组 ， 
第 2 组 为 其 他 种 群 组 ,进行 分 子 方差 分 析 ( 表 8)。 
从 总 群体 分 析 来 看 ,不 同 种 群 间 存在 一 定 的 遗传 分 
化 (Fs =0.131, P<0.001)。 其 中 ,87% 的 遗传 变 
异 来 自 种 群 内 部 。 变 异 等 级 分 析 ( hierarchical 
AMOVA) 分 析 结 果 表 明 ,不 同 组 之 间 存 在 明显 分 化 ， 
占 全 部 遗传 差异 的 25% , 同一 区 域内 不 同 地 理 种 群 





































































































80 


昆虫 学 报 Acta Entomologica Sinica 63 卷 





JL SY2014 
信托 


SY2012 A 
SY2013 


PCoA2 (14%) 





PCoA1 (43%) 





图 2 中 国 东北 地 区 15 个 灰 飞 恒 地 理 种 群 基于 Rsr 值 的 遗传 距离 矩阵 的 主 坐 标 分 析 (PCoA) 
Fig. 2 ”Principal coordinate analysis (PCoA) of 15 geographical populations of Laodelphax striatellus 


in Northeast China based on the genetic distance matrix of F's values 











轴 1 和 轴 2 分别 解释 了 43 % 和 14 % 的 变异 。 Forty-three percent of the variation was explained by axis 1 and 14% by axis 2. 
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图 3 lnP(K) 值 变动 图 和 AK 方法 绘制 的 K 值 变动 图 
Fig. 3 The graph of the value lnP(K) in each K and the number of the inferred clusters (K) by AK method 
1 -15: 分 别 为 15 个 地 理 种 群 Fifteen geographical populations respectively. 
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图 4 基于 9 个 微 卫 星 位 点 的 中 国 东北 地 区 灰 飞 虱 15 个 地 理 种 群 的 STRUCTURE 聚 类 分 析 
Fig. 4 STRUCTURE clustering analysis of 15 geographical populations of Laodelphax striatellus 





in Northeast China based on nine microsatellite loci 
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表 8 中 国 东北 地 区 15 个 灰 飞 虱 地 理 种 群 的 AMOVA 分 析 
Table 8 AMOVA analysis of 15 geographical populations of Laodelphax striatellus in Northeast China 
变异 来 源 自由 度 平方 和 方差 分 量 变异 百分比 P 值 固定 指数 
Source of variation d.f Sum of squares Variance components Percentage of variation P-value Fixation index 
AMOVA 总 群体 分 析 
Global AMOVA analysis 
种 群 间 Among populations 14 63.302 4.522 13 0.001 
种 群 内 Within populations 360 343.600 0.954 87 0.001 0 
分 子 变异 等 级 分 析 
Hierarchical AMOVA analysis 
组 间 Among regions 1 464.471 3;205. 25 0.001 For =0.248 
组 内 种 群 间 Among populations 13 1 094.404 3a1ly 24 0.001 Fs =0.321 
种 群 内 Within populations 360 2 374.293 6.595 51 0.001 Fir =0.489 





间 的 遗传 分 化 占 全 部 遗传 变异 的 24% ,另外 51% 遗 


传 变异 来 自 种 群 内 。 





2.5 遗传 距离 与 地 理 距 离间 的 相关 性 分 析 


如 图 5 所 示 , 对 15 个 灰 飞 乱 地 理 种 群 间 的 遗传 
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距离 与 地 理 距离 进行 相关 性 回归 分 析 , 回归 方程 为 
y=0.0061x +0. 1323(r =0.036, P=0.480), 由 此 
可 见 ,15 个 灰 飞 届 地理 种 群 的 遗传 距离 与 地 理 距 离 
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无 显著 相关 性 。 
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基于 9 个 微 卫星 位 点 的 中 国 东北 地 区 灰 飞 虱 15 个 地 理 种 群 遗 传 距离 与 地 理 距 离间 的 相关 性 分 析 





Fig. 5 ”Regression analysis between genetic distance and geographical distance among 15 populations of Laodelphax striatellus 


in Northeast China based on nine microsatellite loci 





3 讨论 
3.1 种 群 遗 传 多 样 性 

探究 物种 遗传 变异 有 助 于 揭示 物种 起 源 与 进化 
历史 (Lohman et al.,，2008 ) 。 通 常 认 为 ,物种 具有 较 
高 的 遗传 变异 ,说 明 其 可 能 更 具有 较 强 的 适合 环境 
变化 的 能 力 (Eraiksson ef al., 1993 ) 。 等 位 基因 数 和 

















杂 合 度 水 平 是 群体 内 遗传 变异 的 重要 指标 ,同时 也 
是 微 卫星 位 点 遗传 多 样 性 指标 ,但 等 位 基因 数目 容 
易 受 到 样本 量 的 影响 (Maudet et al.，2002 ) ,样本 量 
的 大 小 与 等 位 基因 数 呈 正 相 关 ( 杨 君 等 , 2018 ) 。 等 
位 基因 数 越 多 ,说明 该 位 点 的 遗传 多 样 性 越 高 。 观 
测 杂 合 度 (Ho) 是 一 个 座位 的 杂 合 子 数 除 以 观察 个 
体 总 数 , 它 与 期 望 杂 合 度 (He) 相 比 ,更 易 受 样本 大 
小 等 因素 的 影响 ( 杨 泽 宇 等 ,2007)。Takezaki 和 
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Nei(1996) 提 出 微 卫 星 计算 出 的 期 望 杂 合 度 (He) 
在 0.300 ~ 0. 800 之 间 , 则 可 说 明 群 体 遗 传 多 样 性 


= 
较 高 。 





级 分 析 (hierarchical AMOVA ) 分 析 结 果 表 明 ,不 同 组 
之 间 存 在 明显 分 化 ,同一 区 域内 不 同 地 理 种 群 间 的 
遗传 分 化 , 占 全 部 遗传 变异 的 24% ,另外 ,51% 遗传 























在 本 研究 中 ,15 个 灰 飞 乱 地 理 种 群 375 头 灰 飞 
避 的 SSR 分 析 结 果 发 现 ,9 个 微 卫 星 多 态 性 较 高 , 平 
均 每 个 位 点 对 应 15 个 地 理 种 群 的 等 位 基因 数 介 于 
8 ~45 个 ,平均 为 20 个 ( 表 4)。 说 明 此 9 对 微 卫星 
引物 在 15 个 地 理 种 群 中 所 提供 的 多 态 信息 含量 较 
为 丰富 ,遗传 多 样 性 分 析 结 果 可 靠 。15 个 地 理 种 群 
的 平均 观测 杂 合 度 和 平均 期 望 杂 合 度 分 别 为 0. 548 
和 0.582( 表 5)。 各 地 理 种 群 遗 传 多 样 性 较 高 ,说 
明了 灰 飞 乔 对 环境 的 适应 能 力 较 强 ( 牛 成 伟 等 ， 
2006) 。Sun 等 (2012 ) 开发 了 9 个 灰 飞 乔 的 微 卫 星 
位 点 ,并 对 中 国 3 个 地 理 种 群 (江苏 、 山 东 和 浙江 ) 
进行 微 卫 星 研究 ,结果 显示 这 些 位 点 具有 高 度 多 态 
性 ,每 个 位 点 对 应 的 3 个 种 群 有 13 ~ 30 个 等 位 基 
因 。 蒋 欣 雨 (2014) 开发 了 7 个 微 卫 星 位 点 对 3 个 
灰 飞 乔 种 群 进行 种 群 遗 传 多 样 性 的 分 析 , 得 出 3 个 
种 群 的 观察 杂 合 度 和 预期 杂 合 度 分 别 在 0. 543 ~ 
0.971 和 0.604 ~0.876 之 间 。Sun 等 (2015 ) 对 中 国 
24 个 灰 飞 乱 地 理 种 群 进行 遗传 多 样 性 分 析 的 结 
也 与 本 研究 得 出 结果 相似 ,平均 观测 杂 合 度 Hi 值 
为 0.590 ~ 0.774,uHe 值 为 0.752 ~ 0.810。 总 体 
上 ,我国 东北 地 区 不 同 地 理 种 群 灰 飞 乱 具有 较 高 的 
遗传 多 样 性 。 

3.2 种 群 遗 传 分 化 与 基因 流 

种 群 间 固 定 指数 Fy 反映 了 种 群 间 的 分 化 程 
度 。Fsr 越 接近 0 ,种群 间 越 相似 ; 越 接 近 于 1, 分 化 
越 显 著 。 本 研究 结果 表明 ,不 同 地 理 种 群 的 灰 飞 乔 
具有 一 定 的 遗传 分 化 ,各 地 理 种 群 间 的 遗传 分 化 系 
数 为 0.012 ~0.550 ,各 地 理 种 群 间 的 基因 流 Nm 为 
0.204 ~20.275( 表 6) ,其 中 吉林 种 群 ( 开 ) 与 其 他 
14 个 地 理 种 群 的 基因 流 均 小 于 1 ,遗传 分 化 系数 均 
大 于 0.25 ,说 明 吉 林 种 群 (并) 较 其 他 种 群 有 一 定 的 
遗传 分 化 。 孙 研 涛 (2012) 对 灰 飞 避 种群 的 sj 结果 
显示 总 体 Rs 为 0.0036, 只 有 东北 的 两 个 种 群 有 较 
微弱 的 分 化 ,这 与 本 研究 结果 相似 。 

3.3 种 群 遗传 结构 

通过 构建 UPGMA 系统 发 育 树 、 主 坐标 分 析 
(PCoA) 及 STRUCTURE 聚 类 分 析 进 行 种 群 遗传 结 
构 分 析 ( 图 1 ~4) ,结果 表明 , 灰 飞 乱 具有 一 定 的 谱 
系 遗 传 结构 ,吉林 种 群 ( 柜 ) 相对 其 他 种 群 产生 较为 
明显 的 遗传 分 化 。 此 外 ,AMOVA 分 子 方差 分 析 表 
明 ( 表 8) ,87% 遗传 变异 主要 发 生 在 种 群 内 ,变异 等 






















































































































































































变异 来 自 种 群 内 ,表明 东北 地 区 灰 飞 乱 已 存在 明显 
的 遗传 分 化 ,而 灰 飞 乔 遗 传 分 化 无 一 定 的 地 域 性 差 
异 , 地 理 距 离 对 灰 飞 各 遗传 分 化 没有 显著 的 影响 。 
蒋 欣 雨 (2014) 利用 7 个 微 卫 星 位 点 对 3 个 灰 飞 二 
种 群 进行 种 群 遗 传 多 样 性 的 分 析 , 结 果 显 示 灰 飞 乔 
种 群 分 化 较 小 ,但 东北 黑龙 江 两 个 种 群 与 其 他 种 群 
相 比 分 化 程度 较 高 ,说明 东 北 种 群发 生 了 一 定 程度 
的 遗传 分 化 ,这 可 能 与 灰 飞 异种 群 中 存在 居留 型 导 
致 的 。 这 一 研究 现象 与 本 研究 吻合 。 
东北 夏秋 的 主要 风向 是 西南 风 和 南 风 ,显然 不 

利于 春季 迁 入 的 昆虫 往 南 回迁 ,有 可 能 对 迁 入 东北 
的 昆虫 种 群 具 有 Pied pipe 效应 ( 武 向 文 ,2001 ) 。 
而 国内 报道 灰 飞 蛋 在 东北 地 区 能 够 以 若虫 越冬 ( 林 
志 伟 等 , 2004) 。 利 用 线粒体 基因 (CO7 和 CON) 和 
13 个 SSR 位 点 对 中 国 22 个 灰 飞 恒 地 理 种 群 进 行 种 
群 遗 传 结构 及 种 群 遗 传 多 样 性 的 研究 发 现 ,整个 气 
候 区 的 微 卫星 多 样 性 水 平 大 致 相似 , 由 于 线粒体 基 
因 组 比 核 基因 组 更 易 受 到 影响 ,温带 地 区 的 线粒体 
多 样 性 水 平 明 显 低 于 其 他 气候 区 。 推 测 可 能 由 于 东 
北 地 区 的 水 稻 种 植 面 积 扩大 使 得 灰 飞 乱 大 量 迁 移 导 
致 的 遗传 分 化 (Sun ei al., 2015 ) 。 候 文 术 (2013 ) 比 
较 了 我 国 5 个 地 理 种 群 (郑州 . 瘟 兴 \ 沈 阳 ` 济 末 、 南 
京 ) 中 不 同 灰 飞 乱 种群 的 生态 学 特性 差异 ,从 雌雄 
比 .发育 历 期 有 效 积温 等 方面 发 现 沈 阳 地 区 灰 飞 乔 
种 群 与 其 他 地 区 种 群 有 较 大 不 同 。 从 19 世纪 末 开 
台 , 东 北 地 区 水 稻 种 植 面 积 从 东北 的 东部 地 区 开始 
扩大 , 慢 慢 扩大 到 南部 地 区 和 北部 地 区 ,这 可 能 为 迁 

的 灰 飞 虱 提供 了 一 定 的 食物 来 源 ,因此 东北 地 区 
灰 飞 乱 的 遗传 分 化 可 能 是 种 群 中 迁 飞 型 与 居留 型 的 
分 化 导致 。 
3.4 害虫 防治 启示 与 未 来 工作 

目前 , 灰 飞 乱 广 泛 分 布 于 我 国 各 水 稻 主 产 区 ,该 

种 害虫 的 远 距 离 迁 飞行 为 使 得 我 国 部 分 地 区 灰 飞 乔 
种 群 存在 一 定 的 基因 流 (Sun et al., 2015 ) , 抗 性 基 
因 随 长 距离 迁 飞 而 扩散 ,可 能 导致 杀 虫 剂 抗 性 基因 
的 迅速 传播 。 因 此 ,在 较 大 的 时 空 范围 内 调查 具有 
杀 虫 剂 抗 药 性 的 灰 飞 虱 种 群 ,研究 种 群 中 抗 性 基因 
相对 于 遗传 结构 和 基因 流水 平 的 分 布 ,能 够 确定 其 
对 杀 虫 剂 的 抗 药 性 水 平 。 本 研究 所 使 用 的 灰 飞 恒 样 
品 均 采 目 东 北 地 区 ,下 一 步 仍 需 对 全 国 地 区 的 灰 飞 
虱 种 群 进行 样品 采集 ,针对 全 国 地 区 的 灰 飞 剧种 群 

























































































































































































1 期 姜 ” 姗 等 : 基于 微 卫星 标记 的 中 国 东北 地 


区 灰 飞 乱 遗传 变异 及 种 群 遗 传 结构 分 析 83 





进行 遗传 多 样 性 及 遗传 结构 的 研究 。 
除了 地 理 和 气候 以 外 ,寄主 植物 的 专 化 性 和 杀 

虫 剂 的 施用 也 可 能 影响 灰 飞 必 种 群 遗 传 变异 的 分 布 
( Brévault et al., 2008 ) 。 因 此 ,针对 同一 地 区 灰 飞 乱 
不 同年 际 间 的 种 群 多 样 性 及 遗传 结构 研究 将 有 助 于 
揭示 灰 飞 虱 对 该 地 区 的 适应 性 。 本 文 对 东北 地 区 的 
灰 飞 乔 遗传 多 样 性 及 遗传 结构 进行 了 研究 ,但 由 于 
样品 采集 限制 ,东北 地 区 灰 飞 息 同一 地 区 不 同年 际 
间 的 种 群 多 样 性 及 遗传 结构 尚未 明确 。 此 外 ,近年 
来 以 菌 治虫 的 的 思路 受到 ) 泛 重 视 ， 对 稳 飞 异体 内 

微生物 的 多 样 性 的 研究 ,能 够 为 将 来 以 菌 治虫 打下 
ee et dl.，2016 ) 。 而 针对 灰 飞 乱 在 我 
国 的 种 群 扩 散 路 径 及 种 群 分 化 历史 的 深入 研究 ,对 
灰 飞 乱 的 预测 预报 具有 重要 作用 。 因 此 ,未 来 为 了 
更 好 地 理解 灰 飞 乱 的 种 群 历史 和 种 群 遗 传 结构 及 对 

其 未 来 的 扩散 进行 准确 预测 预报 ,还 需 多 种 手段 进 
行进 一 步 研究 。 
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